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ABSTRACT 

Transient spin-up and spin-down flows inside of 
spacecraft fuel tanks need to be analyzed in order 
to properly design spacecraft control systems. 
Knowledge of the characteristics of angular 
momentum transfer to and from the fuel is used 
to size the de-spin mechanism that places the 
spacecraft in a controllable in-orbit state. 

In previous studies, several analytical models of 
the spin-up process were developed. However, 
none have accurately predicted all of the flow 
dynamics. Several studies have also been 
conducted using Navier-Stokes based methods. 
These approaches have been much more 
successful at simulating the dynamic processes 
in a cylindrical container, but have not addressed 
the issue of momentum transfer. 

In the current study, the spin-up and spin-down 
of a fuel tank filled with gaseous xenon has been 
investigated using a three-dimensional unsteady 
Navier-Stokes code. Primary interests have been 
concentrated on the spin-up/spin-down time 
constants and the initial torque imparted on the 
system. Additional focus was given to the 
relationship between the dominant flow 
dynamics and the trends in momentum transfer. 
Through the simulation of both a cylindrical and 
a spherical tank, it was revealed that the transfer 
of angular momentum is nonlinear at early times 
and tends toward a linear pattern at later times. 
Further investigation suggests that the nonlinear 
spin up is controlled by the turbulent transport of 
momentum, while the linear phase is controlled 
by a Coriolis driven (Ekman) flow along the 
outer wall. These results indicate that the spin- 
up and spin-down processes occur more quickly 


in tanks with curved surfaces than those with 
defined top, bottom, and side walls. The results 
also provide insights for the design of spacecraft 
de-spin mechanisms. 


NOMENCLATURE 

H Angular momentum 

P Pressure 

t Time 

T Temperature 

Re Reynolds number 

v Kinematic viscosity 

p Density 

Q Angular velocity 


INTRODUCTION 

The problem of spin-up from rest is well suited 
to computational analysis because the internal 
flow dynamics have been found to be non-linear, 
and cannot be easily visualized or measured 
experimentally [1]. Hie analysis of spin-up of an 
incompressible fluid in a cylindrical tank has 
been performed in previous studies [1-3]. 
Several of these investigations have successfully 
implemented Navier-Stokes based methods, 
demonstrating significant improvements over the 
results obtain from known analytical methods. 
Hyun, et al. concluded in their work [3] that the 
spin up process is controlled in the following 
manner. Immediately after impulsive start, 
viscous diffusion dominates spin up. Once, the 
Ekman flux is established, the spin up is 
dominated by nonlinear radial advection. 
Eventually, for t/t fmal >0.5, the strength of the 
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nonlinearity diminishes and the linear Coriolis 
force dominates. 

In the current study, the spin-up/spin-down time 
constants and the initial torque imparted on the 
system are of primary interest. Additionally, the 
flow dynamics which control the spin-up are also 
examined. Previous work from Greenspan [2], 
suggests four areas of particular interest. Figure 
1 outlines these areas as: (1) the region in front 
of the shear front, which is at rest, (2) the region 
behind the front, which is in motion, (3) the 
Ekman layer which is a result of Coriolis forces, 
and (4) which is the reverse flow zone that is a 
product of the Ekman layer pumping fluid from 
region 1 to region 2. 


CO 
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Figure 1: General flow pattern during spin up of 
a cylindrical tank: (1) region in front of the shear 
front, (2) region behind the shear front, (3) 
Ekman layer, (4) reverse flow zone 


In the current study, two geometries were 
analyzed; a cylindrical tank, and an oblate 
spheroid tank. Both cases were impulsively 
started from rest and spun up to a constant rate of 
rotation and then spun back down to rest. The 
flow conditions for temperature, pressure, 
rotation rate, etc. were based on estimated 
conditions of a recent spacecraft mission. The 
investigations were performed using a three- 
dimensional unsteady Navier-Stokes code. This 
code is particularly suited for rotating 
components and has been highly successful at 
predicting unsteady phenomena in 
turbomachinery applications. 


NUMERICAL PROCEDURE 


The general code structure is based on a well- 
established compressible, 3-D, unsteady 
turbomachinery flow code [4]. It employs a 
system of overset O-grids and H-grids, with the 
values on the O-H boundaries being updated 
each time step by bilinear interpolation from the 
adjacent grid. The inviscid fluxes are third-order 
spatially accurate and are calculated using Roe’s 
scheme [5]. The viscous fluxes are calculated 
using second-order accurate central differences, 
and the code is second-order accurate in time. A 
modified Baldwin-Lomax turbulence model is 
used for turbulence closure [6]. 

The code contains two options for the 
determination of fluid properties. The first 
option is based on the equations of state, 
thermodynamic departure functions and 
corresponding state principles constructed by 
Oefelein [7]. The second option is based on 
splines generated from the NIST Tables [8], and 
is computationally more efficient. The code 
currently contains fluid property routines (both 
gaseous and liquid) for water, hydrogen, oxygen, 
nitrogen, RP-1, methane, carbon monoxide and 
xenon. 

Computational Grids 

In the current study, algebraically generated H- 
grids were used to discretize the entire flow field. 
The computational grids are shown in Figs. 1 
and 2 for an oblate spheroid tank and a cylinder, 
respectively. Computational grids consist of 
approximately 500,000 grid points with 
clustering near the top and side wall boundaries. 
The average value of y+, the non-dimensional 
distance of the first point off the wall, is less than 
1.0 in both cases. 



Figure 2: Computational grid for an oblate 

spheroid tank 
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Figure 3: Computational grid for an oblate 
cylinder 


Boundary Conditions 

At solid surfaces the relative velocity is set to 
zero, the normal derivative of the pressure is set 
to zero, and the surfaces are assumed to be 
adiabatic. 

A detailed description of the code/algorithm 
development, as well as its application to several 
turbine and pump test cases, is presented in Refs. 
9 and 1 0. 


NUMERICAL SIMULATIONS 


The numerical simulations were performed for 
both geometries using the following flow 
conditions: 


T = 25 C 
P = 1250 psi 
p = 1667 kg/m 3 
v = 5.204 x 1 O' 8 m 2 /s 
Q = 60 rpm 
Re =2.3 x 10 7 

The flow is essentially incompressible, and the 
pressure and temperature (thus the density) vary 
little throughout the spin-up process. The tanks 
were impulsively spun up from rest and the 
subsequent gas spin up was observed. Likewise, 
the rotation of the tanks was impulsively 
stopped, and the tanks were spun down to rest to 
observe the gas dynamics phenomena. 

In Fig. 4, the velocity magnitude contour for the 
cylindrical tank geometry is shown. The shear 
front and reverse flow zones are clearly present. 



Figure 4: Velocity magnitude contours for the 
cylindrical tank geometry 


In Fig. 5, the reverse flow zone is magnified and 
shown with the addition of velocity vectors. 
From this view the Ekman flow along the top 
wall is clearly visible and shown to turn 
downward as it encounters the outer wall and 
develops into a vortex. 



Figure 5: Magnified contours of the Ekman flow 
and the reverse flow zone 

The angular momentum trend for the cylindrical 
tank is plotted in Fig. 6. It can be seen that at 
early times the spin up of angular momentum is 
non-linear. At later times the trend tends toward 
a more linear growth. 
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Figure 6: Angular Momentum vs. Time for the 
Cylindrical Tank Geometry 


To further investigate the driving factors behind 
the trend in angular momentum spin up, the 
turbulent viscosity at the tank equator was 
examined. In Fig. 7, the viscosity factor, which 
represents a ratio of turbulent viscosity to 
kinematic viscosity, is plotted versus radial 
position. At early times the turbulent viscosity is 
very large and dissipates to a much smaller value 
after about 15 seconds. It is also observed that 
the turbulent viscosity is highest near the outer 
wall and rapidly decreases as you move radially 
inward. This indicates that at early times there is 
a strong turbulent transport of momentum 
pushing the shear front inward. As time 
increases, this transport weakens. 



Figure 7: Viscosity factor vs. radial position 

(0=outer wall) at various times 


In Fig. 8, the angular momentum spin-up and 
spin-down trends are plotted for the spheroid 
tank case. The spin up trend in momentum is 
non-linear for t < 0.5 t fina i. Approximately 85% 
of the total momentum transfer occurs during 
this period. Knowing that torque is simply the 
time derivative of angular momentum, it can be 
seen that the initial torque for the spherical tank 
is larger than that for the cylindrical tank. Figure 


8 also shows the spin-down of the tank initiated 
at two different times. In both cases, the general 
trend of spin-down mirrors the shape of the spin- 
up. 



Figure 8: Angular momentum vs. time for 
spheroid tank 


It can be seen in Fig. 9, that the spin-down initial 
torque is slightly larger than the spin-up. Also, it 
is seen that the second spin-down case (red 
curve) has a slightly smaller initial slope than the 
first (blue curve), indicating less turbulent 
viscosity. This is intuitive as the first (blue case) 
begins the de-spin process at time in which the 
flow is in a more unsteady state. 



Figure 9: Comparison of the initial slope of 

spin-up (black) and spin-down (blue and red) 


Figure 10 shows the flow contours in the system 
during spin-up. The zones of reverse flow are 
located at the equator of this geometry and the 
Ekman flow develops along the entire length of 
the outer wall. 
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Figure 10: Velocity magnitude contours for the 
oblate spheroid tank geometry, t=30 s 


A more detailed view of flow near the outer wall 
is shown in Fig. 11. Near the wall, flow is seen 
to be moving toward the equator (down), while 
just away from the wall the flow can be seen to 
be moving away from the equator. The 
development of several vortices is noted, 
indicating that flow behind the shear front is not 
necessarily in solid body rotation. 



Figure 11: Velocity vectors near the outer wall 
of the oblate spheroid tank 

[To Be Added] 

Results for Q = 50 rpm 
Velocity Profiles 
Particle Streamlines 
Boundary layer profile in the 
Ekman layer 


CONCLUSIONS 

A series of numerical simulations have been 
performed to study the unsteady flow 
phenomena and momentum transfer associated 
with the spin-up and spin-down of fuel tanks. 


In the simulation of a cylindrical tank, it was 
found that at early times the angular momentum 
spin-up increases rapidly and in a nonlinear 
fashion. Indications from viscosity profiles 
suggest that during this time the flow is 
dominated by turbulent transport of momentum, 
which pushes the shear front radially inward. At 
later times, the turbulent transport weakens and 
the trend in angular momentum tends toward a 
more linear growth. 

The simulation of the spheroid tank demonstrates 
many of the same dynamics as those seen in the 
cylindrical tank. A significant difference was 
seen in the development of the Ekman layer 
along the entire length of the outer wall, which 
resulted in the reverse flow zones being located 
at the equator of the tank. As a result the 
turbulent transport process is stronger in the 
spheroid tank, resulting is a faster spin-up 
process by approximately 50%. 

Overall, the results of this study suggest that at 
early times the spin-up process is controlled by a 
non-linear transport of momentum. At later 
times the flow is controlled by Coriolis effects 
and the spin-up process tends to be more linear 
in nature. Lastly, the location of the reverse flow 
zones (as controlled by surface curvature) can 
significantly impact the strength of the non- 
linear phase of spin-up. 
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